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Abstract 

The calculation of the dephasing time in electron systems is presented. By 
means of the Keldysh formalism we discuss in a unifying way both weak lo- 
calization and interaction effects in disordered systems. This allows us to 
show how dephasing arises both in the particle-particle channel (weak local- 
ization) and in the particle-hole channel (interaction effect). First we discuss 
dephasing by an external field. Besides reviewing previous work on how an 
external oscillating field suppresses the weak localization correction, we derive 
a new expression for the effect of a field on the interaction correction. We 
find that the latter may be suppressed by a static electric field, in contrast 
to weak localization. We then consider dephasing due to inelastic scattering. 
The ambiguities involved in the definition of the dephasing time are clarified 
by directly comparing the diagrammatic approach with the path-integral ap- 
proach. We show that different dephasing times appear in the particle-particle 
and particle-hole channels. Finally we comment on recent experiments. 

I. INTRODUCTION 

Quantum corrections to the classical Drude formula for the electrical conductivity give 
rise to singular terms in low-dimensional systems [|I]|§. There are two types of such terms. 
The first is known as the weak localization correction (WL) and arises as a result of the 
quantum interference of electron waves in disordered systems. The second is a consequence 
of the enhancement of the electron-electron interaction in a disordered system and is usually 
referred to as electron-electron interaction correction (EEI). In two dimensions, in particular, 
the correction to the conductivity is logarithmic for both the WL and EEI correction. The 
argument of the logarithm contains the ratio of two length scales. The first is the mean 
free path, I, which sets the microscopic scale, beyond which the system behaves diffusively. 
The second length scale differs in the two cases. For WL it is L^, the scale over which 
inelastic scattering starts to destroy the interference effects. For EEI it is usually given by 
the thermal length Ly. In a diffusive system, all length scales correspond to characteristic 
times: the elastic scattering time r = l 2 /D, the dephasing time = L^/D, and the thermal 
time tt = Ti/ksT = L^/D. In general, dephasing occurs either as a result of the interaction 
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of the system with the external environment or as a consequence of the internal electron 
interactions. For example, it is known that an AC-field suppresses the WL. On the other 
hand, while it is clear that inelastic scattering contributes to dephasing, the exact way this 
happens is by far less obvious as witnessed by the recent hot debate in the literature 
Natural questions to ask are: How does the dephasing time that enters the WL correction 
and cuts off the logarithmic singularity depend on inelastic scattering and on the external 
environment? Does dephasing also affect the EEI correction? The two issues need to be 
addressed together if one wants to understand the experimental data in detail. Although 
one can find in the literature discussions of both these issues, there seems to exist no unified 
discussion of them. It is one of the aims of the present paper to fill this gap. 

It was first noted by Schmid || that the inelastic quasiparticle scattering time is en- 
hanced in the presence of disorder. The dephasing time, which controls the WL correction 
was initially assumed to be identical to the inelastic quasiparticle scattering time in a disor- 
dered Fermi-liquid [[/J. According to this analysis, the inverse dephasing time was assumed 
to be oc T In T in 2d, thus predicting a violation of the Fermi-liquid behavior at low temper- 
atures. However, Altshuler, Aronov, and Khmelnitsky || (hereafter referred to as AAK), 
by means of a semi-classical path-integral approach to the WL correction, were able to cal- 
culate the dephasing time directly and predicted an inverse dephasing time proportional to 
T, in contrast to Ref . . Some time later Fukuyama and Abrahams (hereafter indicated 
as FA), re-examined the problem in terms of standard diagrams and calculated the "mass" 
term that develops in the particle-particle propagator. They found the same result as Ref. 
0. The origin of the discrepancy of the AAK and FA results lies in a genuine ambiguity 
in the definition of the dephasing time itself. A review of the attempts to clarify this issue 



made at the time [[TO-flj] may be found in [13]. The issue was further complicated from 



the fact that a dephasing time also appears in the particle- hole channel [JTJJ]. The mass 
term in the particle-hole propagator turned out to coincide with that found by FA in the 
particle-particle propagator. 

More recently inelastic scattering and dephasing have been addressed from a general 
point of view by Stern et al [|TI| and in connection with disordered mesoscopic systems by 
various groups [|16[], in all cases confirming the AAK result. Furthermore, during the last 
couple of years, various experimental groups have observed a saturation of the dephasing 
time at low temperatures |T7]-ET]]. Usually this saturation has been related to heating effects 



or to the presence of magnetic moments. In Ref. JTiJ] , however, these possibilities have been 



excluded experimentally, and the observed saturation seems to be in contrast to the available 
theories. If the observed saturation of the dephasing time is due to an intrinsic mechanism 
of disordered electron systems, this will have dramatic consequences on the localization 
theory, such as preventing an insulating ground state. The above consideration led to a 
fresh re-examination of the theory whose historical development has been sketched above. 
In particular, such an intrinsic mechanism has been claimed to be discovered in Ref. ||. 
However in Refs. this new theoretical result has been questioned and at the same time 
in Ref. || a conventional mechanism has been suggested as a possible explanation of the 
experimentally observed saturation of the dephasing time. The authors of Ref. || concluded 
that there is no need to revise the existing theory of the dephasing time in WL. While we 
completely agree with Ref. || on this latter point, we think that the interpretation of the 
experiments is not yet settled and, especially, the effect of dephasing on the EEI correction 
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is not clear. 

In this paper we discuss dephasing in the particle-hole channel, and in parallel we clarify 
the reason for the discrepancy in the AAK and FA calculations. The difference between the 
results may be traced back to a different definition of the dephasing time. In the FA case, 
the dephasing time is assumed to be the mass of the particle-particle propagator as due to 
the self-energy corrections. In this type of corrections, one does not include diagrams which 
connect the upper and the lower electron lines entering the particle-particle propagator. In 
the AAK approach these diagrams, i.e. the vertex corrections, are taken into account. The 
inclusion of these vertex diagrams in the diagrammatic approach of FA leads to the AAK 



result ||22|| . In the particle-hole channel, on the other hand, we find that vertex corrections 
are negligible and the dephasing is entirely determined by the self-energy. 

To see this we will adopt a real time-formulation based on the Keldysh technique which 
allows us to discuss WL and EEI corrections in a unified way. Furthermore, in the real-time 
formulation non-linear effects, for example those due to an external field, can be incorporated 
more easily. In order to address the issue of how dephasing arises in the particle-particle and 
in the particle-hole channels, it is instructive to start considering the situation of an external 
field. This analysis, besides providing a relatively simple physical situation to analyse has 
also a genuine experimental interest. For example, there exists no satisfactory theoretical 
explanation for the non-linear field effects observed in Refs. p3"-27fl. To this end we will 
devote a substantial part of our analysis to the nonlinear effects on WL and EEI corrections 
for which we shall find an unexpected difference in their dependence on a static electric field. 

The plan of the paper is the following. In the next section, by means of the Keldysh 
formalism, we derive an expression for the additional current in a disordered conductor due 
to both the WL and EEI corrections. The formulae obtained are valid in the presence of 
an arbitrary time-dependent external field. This result allows us, in section III, to discuss 
the WL and EEI corrections in the presence of both a DC and an AC electric field. We 
show that, in the case of a DC field, the EEI correction is suppressed, in contrast to what 
happens for the case of WL correction. An AC field, on the other hand, suppresses both 
the WL and EEI correction, but with a different dephasing time. In section IV, we will 
review the FA calculation of the inelastic scattering time in two dimensions and show how 
the AAK result may be obtained from it. In section V the path integral approach to the 
dephasing time of AAK is briefly outlined and connection is made with the diagrammatic 
analysis. Special attention is paid to how the cancellation of the infrared singularity occurs. 
We extend this approach to a calculation of the dephasing time in the particle-hole channel. 
Finally, section VI is devoted to some discussion of the experimental results. In Appendix A 
we give technical details concerning the calculation of the dephasing time in the particle-hole 
channel. 



II. QUANTUM CORRECTIONS TO THE CONDUCTIVITY 

In a weakly disordered metal, quantum corrections lead at low temperature to deviations 
from the Drude-Boltzmann theory of transport. The weak localization (WL) corrections are 
due to electrons diffusing along closed paths, where quantum interference causes an enhanced 
back-scattering probability. The correction to the current in the sample is given by (T|] 
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AJwlW = -e 2 Dr- f° dijCj^f (r, r)E(t - V ), (1) 

where D is the diffusion constant, r the elastic scattering time and C* is the so-called 
cooperon or particle-particle propagator. Its graphical definition is given in FigJI]. In the 
presence of a vector potential the cooperon is given by the solution of the differential equation 

{2j; + D [-iV + eA(r, t - § ) + eA(r, t + § )] 2 } C^, (r, r') 

= ±<%-7/)5(r-r'). (2) 

Here and below we use units where % = ks = 1. In addition to WL there are quantum 
corrections to the conductivity due to the interplay of disorder and electron-electron inter- 
action (EEI). In the long- wavelength, low- frequency limit the electron motion is diffusive, 
leading to an enhancement of the effective interaction, V(q, u) — > V(q, \uj + Dq 2 ) 2 , 

and results in either an enhanced or suppressed conductivity in the case of the direct and 
exchange contribution of the interaction, respectively. 

In the following we derive these corrections to the conductivity in the presence of an 
external electromagnetic field. Since one has to calculate a time-dependent, non-linear 



response function, we work within the Keldysh formalism, using the notation of Ref. |28 
This has the advantage that no analytic continuation is necessary at T 7^ in order to get 
the physical response functions. In particular the two-particle propagators like the cooperon 
are directly defined in terms of a retarded and an advanced Green's function, and there is 
no restriction on the energies of these Green's function in contrast to what is found within 
the Matsubara formalism. 

The starting point is the general equation for the current in the presence of both inter- 
actions and an external field. For fixed disorder configuration, the current in ^-dimensions 
is found from the Keldysh component of the electron Green's function, 

A jEEiW=i/^^A^(tk;tk). (3) 

In order to simplify the notation, we only write down the time dependence of the Green's 
function in the following. To first order in the interaction the correction to the Green's 
function is given by 

AG K (t,t) = G R {tM)^ R {tiMG K {t2,t) 
+G K (t,t 1 )Z A (t 1 ,t 2 )G A (t 2 ,t) 

+G R (t,t 1 )Z K (t 1 ,t 2 )G A (t 2 ,t), (4) 

where G is the non-interacting Green's function in the presence of both disorder and an 
external field, S is the self energy, and one has to integrate over t\ and t 2 . It turns out, 
that only the first two terms contribute to the quantum correction to the conductivity, so 
the term G R T> K G A will be neglected from now on. The first two terms in eq.(H) correspond, 



upon impurity averaging, to the diagrams of Fig.2 of the first of Refs. ||29|| . To lowest order 
in the interaction the retarded self-energy is given by 
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(h,t 2 ) = ^(G K (t 1 ,t 2 )V R (t l ,t 2 ) 
+G R (t 1 ,t 2 )V K {t 1 ,t 2 )) 



(5) 



where V R and V K denote the retarded and Keldysh component of the interaction. Singu- 
lar corrections due to the interplay of disorder and interaction arise, after averaging over 
disorder, from singular vertex corrections, i.e. diffusons, which appear on the vertices that 
connect an advanced and a retarded Green's function. Corrections to AG K where both 
vertices are singular are much larger than the contributions with only one singular vertex. 
As a consequence, these various, less singular, terms are neglected. In particular, the third 
term in eq.(f|), and the second term in eq.(|5]) give rise to only one singular vertex and are 
therefore neglected. 

The structure of the impurity averaged Keldysh Green's function is then 



AG K (t,t) = -G H {tA)UtiA-M)G*{tlt'i) 



*V H (ti - t 2 )A 2 (4 £ 2 '; h)G«(t' 2 , t) - c.c 



(6) 



where Go is the disorder average of G and we have indicated by A x and A 2 the vertices 
dressed by disorder. Notice that AG K is purely imaginary, and the current in eq.(Q) is 
real. In order to make the structure of this equation clear, we represent AG K (t,t) and the 
vertices Ai, A 2 graphically in Fig.[2[ An external electromagnetic field affects all the Green's 



functions and all the vertex functions in this expression. However, it is known |29] that the 
sum of the diagrams with two diffusons, i.e. the sum of the terms with an electromagnetic 
field insertion in the Green's functions, cancel. Diagrams which contribute to the quantum 
corrections to the conductivity must have at least three diffusons. In our formalism these 
diagrams are generated by inserting the electromagnetic field in the vertex corrections Ai 
and A 2 . The singular part of the vertices is given by 



A l (t[,t'(;t l ) 
A 2 (£ 2 ; ^ 2 ; h) 



(G fl (fi,*i)G A (ti,t?)> 



2nN r 
1 

27rN r 



{G A (t' 2 \t 2 )G K (t 2} t' 2 )) 



(7) 
(8) 



where the brackets denote the impurity average. N Q is the single-particle density of states. 
A 2 can be easily expressed in terms of the diffuson, 



(9) 



with 



+ D 



-iV + eA(r,t- 



2' 



eAfr, t 



Hit - t')5(r - r>) 



(10) 



The graphical definition of the diffuson is given in Fig.|l|. In comparing the diffuson and the 
cooperon in Fig.|I|, it is worth recalling that the time evolution of the diffuson is associated 
with the center-of-mass time of the particle-hole pair, whereas the relative time controls 
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the evolution of the particle-particle pair of the cooperon. It is also useful to compare the 
different way an external vector potential enters in the equation for the diffuson and the 
cooperon. In the case of the diffuson, the vector potential is felt only if the particle and 
hole are delayed with respect to one another, i.e. when r\ ^ 0. In this respect, processes, 
which are characterized by the simultaneous creation of a particle-hole pair in a given point 
in space, are not affected by an external vector potential. This is the case for the density 
correlation function which is relevant for screening the Coulomb interaction and justifies 
the fact that the screened Coulomb interaction, in eq.(15) below, does not depend on the 
external vector potential. Also the particle-hole pair, which appears in the vertex A 2 , defined 
above, is composed of particles created simultaneously. A 2 therefore does not depend on the 
external vector potential. By contrast the singular part of the vertex Ai also involves the 
Keldysh component of the disorder averaged Green's function, when it is expressed in terms 
of a diffusion propagator see Fig.|l]c. The Keldysh function decays only slowly in time and 
introduces a time delay between the creation of the particle and hole. We approximate the 
Keldysh function by its equilibrium value 

G${u) = tanh(cj/2T) [G r {uj) - G*{uj)] , (11) 

where T is the electron temperature. Notice that in strong electric fields where the electron 
temperature increases due to heating, the electron distribution function may differ from the 
Fermi function, and eq.([TI|) may break down. Such effects are neglected in the present work. 
Going back to the real-time formalism, we Fourier transform G$ . The Fourier transform of 
tanh(c<j/2T) is S(t) — iT/ sinh(7rTt), which has to be convoluted with the advanced Green's 
function. The delta function S(t) can be neglected since it eventually leads to a term 
proportional to 6(t 2 — ti)V R (ti,t2) in the self-energy S(ti,t 2 ) which is zero. Ai is then 
found as 

M*i,#*i) = /d^^ (12) 

and definitely depends on the external field due to the time delay rj in the particle-hole 
propagator. Calculating the current requires also momentum integrations. If the interaction 
V transfers a momentum q with qvFT <C 1, the k-integration over the Green's functions 
gives 

^*<tf(t - 4 k)<#& - V ~ 4 k - q)Gf & - t, k) 

= An N t 2 J eD^t - - rj - 1' 2 ). (13) 

smh(7Ti rj) 

The temperature dependent factor is due to G K (t' 2 — t), with t' 2 — t = —rj. Note that the 
vertex corrections Ai and A 2 do not depend explicitly on k but only on the difference of 
momenta k — (k — q). Finally, we obtain the correction to the current due to the interplay 
of interaction and disorder as 

AN r 2 e f , , , ^ _ / ttT \ 2 



/ n 4iV r e r , , , ^ _ 
AjEEi(t) = / dtidtadr/^Dq 

7T J 



sinh(7rT?7) 

t-T)/2,ti-ri/2^ 



V R ( tl - t 2 , q) At^(q) A%/2,W 2 (q)- (14) 
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Here for simplicity we assumed a homogeneous external field. In an inhomogeneous field the 
diffuson is a function of two momenta, and the current is determined from -D(q, q). In the 
case of long range Coulomb interactions, the expression for the current simplifies. In the 
relevant low momentum, low frequency region the screened Coulomb interaction is given by 

p . , 1 — iu + Dq 2 

The frequency dependent numerator cancels when multiplied with the diffuson, so that 
the product is frequency independent, and transforms into a delta function after Fourier 
transformation. Thus the convolution of the interaction with the first of the two diffusons 
appearing in eq.(|14D gives 

J dt 2 v R (h - 1 2 , q ) At,7°,(q) = ^r^(*i " 1 + v)> ( 16 ) 

and the expression for the current becomes 

2re ^ f°° . q / nT 



AjEEl(t) 



sinh (71^77) J 



This equation will be the starting point for our further investigations. We cut off the time 
integration with r since the diffusive regime requires that all times are longer than r. We 
see that as with Ajwl, Ajeei is also related to the propagation of two particles. However, 
whereas WL is related to the propagation around a closed loop, i.e. the probability of return 
C(r, r), this is not the case for EEI, due to the prefactor q/g 2 in the momentum integration. 
Also the time dependencies are different. In the case of WL two particles start diffusing at 
t — 7] at rj and return at t at the same point. From the time dependence of the diffuson 
in eq. (|PTD , we conclude that, in the case of EEI, the first particle starts at t — 2r/ at and 
arrives at t — r\ at Yf. The second particle starts at t — rj at and arrives at t at the final 
point rj. 

Equation ( |P7D allows us to obtain the results known in the literature. For instance, the 
interaction correction to the conductivity is found by expanding eq.(|I7]) to first order in the 
electric field. By using eq . ([10|) one finds the diffuson to linear order in the electric field as 

^^(q) = \z- Dq2 " (l + 2/^eq • E V 2 + • • •) (18) 

and one arrives at 

4e 2 Dl r , ( *T V V f 1 \ d/2 

where d is the dimension. The hyperbolic sine provides a cutoff for large times rj. By 
approximating the temperature dependent factor with 1 for r\T < 1 and for rjT > 1, 
one finds immediately the temperature dependence of the correction to the conductivity 
as Acteei oc 1/VT, InT and y/T in d = 1,2,3. In one and three dimensions, the correct 
prefactor is only found by integrating eq . (|19D with the result 
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f -2.5 - 2 4jDjf (d=l) 
Aa EEI ~ ln(l/Tr) (d = 2) (20) 

{ 1.8. £s^T/D (d=3), 

where we subtracted temperature independent terms. Here we only considered the contribu- 
tion to Aueei from the exchange diagram shown in Fig.Q. A correction to the conductivity 
with the same temperature dependence as in eqs. (|20|) is also found from the direct (Hartree) 
diagrams, as shown in the literature These will then change the total amplitude of the 
quantum correction to the conductivity due to the interplay of disorder and interaction. We 
will not consider this effect here and below. 

For completeness, we recall also the results for the WL correction to the conductivity in 
a weak static field. The cooperon is given by 

1 / 1 \ d/2 

Here a phenomenological dephasing time has been introduced in order to cut off the 
infrared singularity. The microscopic origin of such a dephasing time will be discussed later. 
By inserting eq.([H|) in eq.([|), one arrives at the well known expressions 

Aa WL = -jgrlnfo/r) (d = 2) (22) 
{^Jl/D^ (d=3). 



III. DEPHASING BY A HOMOGENEOUS ELECTRIC FIELD 

The quantum corrections to the conductivity may be suppressed due to various mecha- 
nisms. We start by studying the effect of a homogeneous external electric field. The effect of 
an electromagnetic field on WL has been described in the literature Here we consider 
in addition the effect of such a field on EEL We will first consider the non-linear response 
to a strong DC electric field, then we will discuss the response to a DC field in the presence 
of a microwave field. 

From eqs. and ( |TTD both the linear and the non-linear response due to quantum inter- 
ference can be determined. Note that heating effects are neglected in the present discussion. 
In a static electric field, the vector potential is given by A(i) = — Et. The cooperon is then 
found as 

-£>(q- 2eEt)V 

•/ 

(23) 

i.e. the electric field drops out. As a consequence, weak localization is not affected by a 
static electric field HQ . 



Cj _„(!-, r) = — 5>xp 



2_, exp —Dq 2 i] 



S 



There are, however, non-linear terms in the interaction correction to the current. The 
diffuson with the time arguments as needed for the calculation of the current is given by 

A%/ 2lt -3/2,(q) = \ expt-^q - eE V ) 2 V ). (24) 

In contrast to weak localization, the interaction correction is not related to a simple sum 
over the momentum entering the two-particle propagator, since the factor q/g 2 has to be 
taken into account. As a consequence, by shifting the q integration by e~Er] in eq.(p4[) it does 
not simplify the calculation of the current, and especially there are non-linear electric field 
effects. One can easily estimate beyond which electric field the non-linear effects become 
sizeable. Equation ( p4|) can be expanded in the electric field under the condition eEr] < q. 
On the other hand, as far as the interaction correction to the current is concerned, only 
times i] < 1/T and momenta q > 1/Lt = (T/D) 1 ^ 2 contribute effectively. The electric 
field can thus be considered small when eE < T/Lt which is equivalent to the condition 
D(eE) 2 = Tq < T 3 . Non-linear effects are expected to dominate for eE > T/Lt, i.e. when 
the voltage drop over a thermal length exceeds the temperature. For the sake of definiteness, 
let us consider eq.([H]) in two dimensions. Performing the q-integration, the current reads 

e 2 f°° drt ( nTri \ 2 sinh(T n V/2) , , , 

This does not lead to a simple linear dependence of the current on the electric field and 
thus non-linear field effects are definitely present. In the high temperature limit T ^> To, or 
weak electric field limit, it is the temperature dependent factor that cuts off the integral at 
large time rj, whereas the T dependent factor varies only weakly and is approximately one. 
Including the first non-linear correction, the current is determined as 

Aj EEI - -^E (ln(l/Tr) - 1-62^^) , (26) 

where 1.62 is the approximate numerical factor from the %]- integration. 

In the limit of large electric field we approximate the To dependent factor as 

sinh(T V/2) f 1 forTo^l 

T V/2 Pl oV 1 ' ~ \ 1/(W) for^ > 1 V ( > 

and then we split the 77-integration in eq.(]T7|) as 



T 



r] A/T TqT] 4 \smh(irTri) 



The second integral gives only a small contribution to eq.(p8|) provided that T r 1, and 
the value of the current at large fields is found to be 



e 2 



Aj EEI = -E— ln(l/T r). (29) 
Considering now the current as a function of temperature, one will observe the InT be- 



havior for T > T while the current will saturate to eq. (|29|) for T < T . For illustration, we 
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plot in Fig.|3] the non-linear conductivity as a function of temperature obtained by numeri- 
cally integrating eq . (p5|) . The saturation of the conductivity at low temperature is apparent. 
In addition Fig.|3| also shows the conductivity in the presence of a slowly (in time) oscillat- 
ing field Eac which is superimposed to the current driving field E. This situation will be 
examined more closely below. In Fig.|3] the conductivity is determined from the average of 
the current over one period of the oscillating field divided by the static field. The strength 
of the oscillating field is Eac — 10-B. At high temperature, the current is not affected by the 
oscillating field. Since the oscillating field is much stronger than the DC field, the current 
saturates at much higher temperature than in absence of the AC field. In more than one 
dimension one has to distinguish the two cases E _L Eac an d E 1 1 Eac- As it is apparent from 
the figure, a parallel AC field saturates the current at higher temperature. If the frequency 
of the AC field is not too large (uT < 1), the effect can be understood from eq. (|29|) by 
considering the current in a strong electric field E total = E + Eac(^) and expanding in the 
static field 



e 2 



Aj E EI = t 

h 



D7T 07T -&AC 



(30) 



Upon time averaging the contribution, in the first term, not proportional to the DC field, 
vanishes. The difference between parallel and perpendicular AC and DC field stems from 
the last term and leads to Aa\\ — Aa± = (2/37r)(e 2 //i) in agreement with the numerical 
result in Fig.|3| in the low temperature region. 

We find analogous results in one and three dimensions. Figure |] depicts the current 
as a function of temperature obtained by numerically integrating eq. flTT| ) in one dimension. 
Again we studied the current in the presence of a static field only, and the time averaged 
current in the presence of both a static and an oscillating field with Eac = 10J5. As in 
two dimensions, one observes a saturation of the current in the low temperature region. 
The temperature scale, where the saturation occurs is roughly the same in one and two 
dimensions. 

We consider now specifically the response in the presence of time-dependent electric 
fields. The diffuson is then given by 

Dt-T]/2,t-3r)/2 (l) = r 

x exp {- J° v dhD (q + eA(t + h) - eA(t + t x - v)f} • (31) 

We start with fields which are only slowly time dependent and discuss high frequency fields 
later. As we mentioned before, only small time delays, i.e. r] < 1/T, are relevant for the 
calculation of the interaction correction to the current. Therefore, when the vector potential 
varies only slowly on this time scale, it may be expanded in the equation above, 

A(t + ti) - A(t + ti -77) « -77E(*), (32) 

and one finds that the current instantly follows the electric field in this case. The conduc- 
tivity in the presence of a slowly oscillating field in Figs.|3| and ^ was determined using this 
approximation. 

In the case of a fast time dependence of the vector potential, the approximation ([32]) is not 
sufficient and the full time dependence has to be taken into account. For A = Eac cos(u;i)/a; 
the relevant combination is 
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A(*i + r]/2) - A(ti - r]/2) = -2E AC Bin(w*i) sin(^/2) /w. (33) 



In the high frequency limit the sine-functions are of order one, so that compared to (p2|) the 
vector potential is reduced by a factor u)r\. The linear response to a high frequency electric 
field, for instance, becomes temperature independent for u 3> T, since the main contribution 
to Acteei, eq.(P7|), comes from 77 < 1/uj <C 1/T and processes with ur] > 1 are practically 
cut off. The conductivity is then found as Acteei oc uj 1 ^ 2 , ln(u;r), uj~ 1 / 2 in three, two and 
one dimension, as it is discussed in the literature Q. 

We now concentrate on the effect of a high frequency radiation on the response to a static 
field. In Fig.[|we plot the static conductivity obtained by integrating eq . fll7D in the presence 
of both a static and a high frequency field in two dimensions. The static and AC field are 
parallel. The temperature is T = lOTo with T 3 = D(eE) 2 . We find that the AC field reduces 
the quantum correction to the conductivity. The effect is strongest for low frequencies, 
uj < T, whereas in the limit of high frequencies the effect becomes weaker. In the relevant 
case of weak field and low frequencies, the correction to the static conductivity can be 
determined from ([|| with the result that a EE i(^Ac)-^EEi(0) = 0.078-(e 2 /2n 2 )D(eE AC ) 2 /T 3 
if the AC and DC field are parallel and one third of this for perpendicular AC and DC fields. 
Concerning the high frequency limit we argued after ( |33"D that the strength of the electric 
field in a high frequency is reduced by a factor rju which corresponds to uo/T if rj ~ 1/T. 
We have verified this estimate numerically and analytically and found that for week field 
the correction to the conductivity changes into cr EE i(E AC ) — °eei(0) oc D(eE AC ) 2 / (u 2 T) in 
the high frequency limit. 

To conclude this Section, for the sake of completeness and comparison, we recall the 
results for the weak localization correction in a microwave field [|TJ. Notice that in the 
cooperon we need the sum of the vector potential at times t ± rj/2, 

A AC (ti + V/2) + A AC (ti - ri/2) = 2A AC cosM) cos(cur//2), (34) 

instead of the difference that enters the diffuson. We average the cooperon over one period 
of the AC field and sum over momentum. After some algebra one arrives at 



2e 2 D f°° dr? „ u 
WL vr (AnDy/ 2 Jr V d / 2 




u 3 



I. 



UJ 3 



(35) 



sin(2x) _ 2 sin^\ 
2x x 2 J 

where Tq = D(eE 'ac) 2 and Iq(x) is the imaginary argument Bessel function. 

The external microwave field does not affect the cooperon for small rj, but does so for 
large rj. The estimate of the typical time scale r\ = r A c, at which the external field modifies 
the ^-integration in the equation above, is found from the condition that the argument of 
the exponential in the second line of eq.(B5j) is of order one, and one gets 



f 2D(eE AC ) 2 /u 2 for D(eE AC ) 2 « u 3 

/TAC \ [D(eE AC ) 2 uj 2 /180] 1/5 for D(eE AC ) 2 » u 3 . [ } 
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Strong non-linear effects of the microwave field associated with WL will be seen when tac 
is smaller than the dephasing time r^. For a fixed strength of the microwave -Eac> the time 
scale tag becomes large when lo — > 0, tag oc uj~ 2 ^ 5 , which agrees with the observation that 
a static field does not lead to nonlinear effects in WL. 



IV. THE INELASTIC SCATTERING RATE 

In this section we calculate the "mass" in the cooperon and diffuson propagator. This 
section does not contain new material, since we will reproduce the results of Refs. |), 10 , 14 



However we think that this section contains useful information in order to understand the 
different approaches to the calculation of the dephasing time. 

In the presence of time reversal invariance, the particle-particle propagator may be ob- 
tained from the particle-hole propagator by reversing the direction of one of the electron 
lines. As a consequence, the mass term that develops in the particle-particle propagator 
also enters the particle-hole propagator. A mass term in the particle-hole channel could 
seem to imply, at first sight, a violation of the particle conservation. However, this is not 
the case because such a mass term eventually disappears in the physical response function, 



as it has been shown by the detailed diagrammatic analysis of Castellani et al. [fL4|| . We 



now review briefly the derivation of how such a mass term arises. For our convenience, we 



adopt the notation of Ref . |14j , and restrict the discussion to the two-dimensional case. The 
two-particle propagator is denoted by L(Q, fl) and may represent a cooperon or diffuson. 

According to the standard diagrammatic technique, a mass term for a propagator may 
be obtained by considering the appropriate self-energy. The self-energy diagrams for the 
two-particle propagator are those labeled (a)-(d) in Fig.|| By denoting with E(Q, fl m ) the 
contribution of these diagrams, one can use the Dyson equation to re-sum their infinite series 
and obtain the two-particle propagator L(Q, fl) as, 

L(Q ' Qm) = sw^wj- (38) 

Successively, by expanding the self-energy in powers of fl and Q 2 , one may identify the 
quantum corrections to the diffusion constant, the density of states, and the frequency. This 
has been done explicitly in Ref. P0[ . The mass term, which will be called l/r ine i, is found 



from the self energy at zero momentum and frequency. The only part of the graphs in Fig|| 
which is not manifestly proportional to Q or fl, is from graph (d) and is given by 

— = -2T J2 £nq^)£(q^ + ^n)- (39) 

Tinel e n <LO v <e n +Q m q 

Here, e n = 7rT(2n + 1), lo v = 2vitT, and fl m = 2mirT are Matsubara frequencies with 
e„, < and e n + fl m > 0. fl m is the external frequency, which after analytical continuation 
ifl m — > fl + i0 + , is sent to zero. In this limit, the number of addenda in the Matsubara 
sum of eq.(^) vanishes. However, the functions in the sum contain a branch cut after the 
analytic continuation and eq.(|59"D yields a finite result in the fl — > limit. 

After the analytic continuation with ifl m — > Q + i0 + and ie n — > e + i0~, by taking the 
limit Q, e —>■ 0, (i.e., particles at the Fermi surface), and observing that the sum of the Fermi 
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and Bose function f(u) + b(u) = [coth(u;/2T) — tanh(u;/2T)]/2 = 1/ sinh(u;/T) is an odd 
function of u>, one gets 

1 2 r°° 1 

— = — duj ■ w m £ ReL(q ' -Mimvfa, ( 4 °) 

Considering the screened Coulomb interaction as in eq . fll5|) , Iml^q, — iu) = —uj/(2N Dq 2 ), 
the integral above is well behaved in the ultraviolet limit both in frequency and momentum. 
Furthermore, due to the presence of the hyperbolic sine, the main contribution to the inelas- 
tic scattering time comes from frequencies u> < T, physically corresponding to the fact that 
there are no available excitations with larger energy. As usual, at a given temperature, only 
real excitations (u < T) generate lifetimes, while virtual excitations (u > T) are immaterial. 
On the other hand, the integral fl40|) is infrared divergent in one and two dimensions. To 
cure this divergence we insert the inelastic time self-consistently in the propagator L(q, u), 
and also calculate the screened Coulomb interaction more accurately. In two dimensions one 
finds 

V(q, = —— — — , 41 

2N Q q Dqn — iu 

with k = 4ire 2 N , the inverse screening length. Since the screening length is usually much 
shorter than the mean free path, eq.(|4l|) and eq . ( 15 ) coincide in almost all the relevant 



parameter range, but not in the extreme low momentum region, where there is no singularity 
in ImV^(q, to) according to eq . (|4l|) . By exploiting the fact that the main contribution to the 
frequency integral in ( [4"0"D comes from the region |u;| <T, we may write in two dimensions 

1 T f T DK 2 q 



I dq duj- 

JO J-T L 



T ine i 2lT 2 N Jo J -T LO 2 + (Dhiq) 2 

x D( l 2 + V^g (42) 

[io 2 + (Dq 2 + l/r inel ) 2 Y 1 ' 

The frequency dependence of the integrand is controlled by the two Lorentzian factors which 
have typical scales 1/T in ei and Duq. Depending whether Duq > 1/T in ci or Duq < 1/Ti nc i, one 
can safely neglect the frequency dependence of the first or second Lorentzian. Performing 
then the uo integration, we find 



T 



( dq— arctan[T/ (Dq 2 + 1/ T inc i)] 

Ja>an L>Q 



Tinel 7T 2 A^o Jq>Q0 Dq 

T rio 

+ — — / dg/tTineiarctan(T/D/tg) (43) 

TT Z N JO 

where go is given by Dq n = 1/T ine i- For 1/T ine i < T, the dominant contribution to the 
inelastic scattering rate comes from the momentum region qo < q < (T/Z)) 1 / 2 , and we find 

' ^ MD^Trl^^-^—logiD^DN.f/T] (44) 



T^ei 4nDN y meu AnDN 
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This is the FA result for the inelastic scattering rate we mentioned in the introduction. The 
AAK result for the dephasing time in the weak localization can be obtained within the same 
calculation, but with a different cut-off procedure of the low q singularity. One may argue 
that the Coulomb interaction cannot contribute to dephasing on distances which are larger 
than the phase coherence length itself, and therefore use Dq 2 = l/r^ as a low momentum 
cut-off in order to determine the phase coherence length. One then finds 

Tr^ HTT * )x ^m- ^ DN ^ (45) 

There remains the task to justify on a more formal level, the origin of the above cut-off 
procedure. Fukuyama suggested |L0| to calculate S(Q, Q) at finite Q, so that the singularities 
in the propagator L(q, to) and the interaction ^(q, to) in the calculation of the scattering rate 
appear at different momenta, and their product is less divergent. However this procedure 
does not cure the 1/q 2 divergence of the interaction and thus cannot justify the difference 
between the FA and AAK results. In the next section we will demonstrate how the low 
momentum cut-off arises in case of weak localization within the path-integral approach of 
AAK [||. We will also see that in the case of the interaction effect, no such cancellation 
occurs, and the dephasing rates are different in the two cases. 



V. DEPHASING DUE TO COULOMB INTERACTION 

In the previous section we calculated the inelastic scattering rate due to the Coulomb 
interaction. We showed that the "self-energy" E(Q, u), which is given by the sum of the 
diagrams shown in Fig.|6], is responsible of various renormalization effects associated to the 
Coulomb interaction. In particular, inelastic scattering arises from only one type of these 
diagrams in the low energy region (\uj\ < T) and within the Matsubara formalism is re- 
lated to the existence of a branch cut contribution. In the present section we address the 
determination of the dephasing time, the time scale over which the quantum correction to 
the conductivity decays. In order to do this calculation one needs to include also vertex 
corrections, i.e diagrams of the type shown in Fig.|6|e. 

The analysis is strongly simplified if we restrict ourselves to frequencies |u;| < T from 
the beginning. In this limit one neglects processes related to the diagrams a,b,c and a part 
of d in Fig.|| In the case of the inelastic scattering time, as shown in the previous section, 
this is immaterial since it only takes contributions from \u\ < T. Hence, by restricting the 
energy transfers to frequencies lower than the temperature some interaction effects are lost, 
but those relevant for inelastic scattering and dephasing are taken into account. 

It is useful at this point to switch again to the real-time formalism. In the limit of 
low frequencies, the Keldysh component of the electron-electron interaction given by V K = 
coth.(u/2T)(V R — V A ) becomes dominant due to the coth(co>/2T) and the retarded and 
advanced parts can be neglected. The Keldysh component of the interaction has the same 
structure in Keldysh space as an external field, and thus the interaction can be treated 
formally as a fluctuating external field. At the end the fluctuation dissipation theorem is used 
to relate the field fluctuations to the internal electromagnetic field at thermal equilibrium. 

With this in mind we start with the calculation of the dephasing rate for WL within 



a path-integral approach, following the literature |q,pj|. The method allows us to include 
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the vertex corrections in an elegant way. Then it will be straightforward to calculate the 
dephasing rate for EEI. 

The particle-particle propagator obeys a diffusion equation whose solution can be written 
in terms of a path integral 



^„,(r,r') = -^- / ' Vv tl exp-(S + Si) (JO) 

ZT Jr ,=r' 

v 



with 



s " = /, *p < 47) 

Si = -ie /"W tl ■ At(r tl ,ti). (48) 

J 7}' 

The vector potential A t (r, t-y) = A(r, t + h/2) + A(r,i — ti/2) is defined as the sum of 
the vector potentials seen by the two particles entering the propagator. To make contact 
with the diagrammatic approach, the potential at times t±ti/2 corresponds to the retarded 
(advanced) electron lines in the diagram. The vector potential describes the effect of the 
electromagnetic field generated by the other electrons. In the case of a fluctuating field with 
only quadratic correlations, one can average over the field fluctuations, leading to 

Si —> S- m t = 

£ JJ dfi JJ dt 2 r t \f£ (Al(r tl ,ti)Ai(r t21 t 2 )}. (49) 
It is useful to define the average over diffusive paths as 

([•••])«, = j Pr t [-..]exp(-^ )/ J Dr t exp(-5o), (50) 

which allows us to write the cooperon in the presence of the fluctuation field as the product 
of the unperturbed cooperon times a perturbation which has to be averaged over diffusive 
paths: 

C^,(r, r') = i- fj^, Vr h exp(-S )(exp(-S int )) d . (51) 



Following Ref. |TT| we then make the approximation 

(exp(-S , int ))d -> exp(-(S , int )d). (52) 

In the diagrammatic language this approximation is equivalent to selecting a certain subset of 
graphs. As realized in Ref. jlTJ one selects diagrams where interaction lines do not intersect 
and thus correspond to graphs (d) and (e) in Fig.[| In some cases (Si n t)d is linear in r\ — rj' 
and then the dephasing time may be easily read off, (S^d — (rj — r/')/2r^. In general we 
determine the dephasing time from the condition (S- mt )d — 1 for rj — —rf = t$. We may 
then write 
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5, 



int 



"T -"7 q, w 

xexp (iq- (r tl - r t J) 



x2 



cosicj — - — ) + cosfu; — - — 



(53) 



where (i'A J ) qiIJ is the electromagnetic field correlator in reciprocal space. The cosine term 
in t\—t2 refers to correlation of vector potentials seen by the same electron at different times, 
while the term in t% + t 2 represents the correlation of the vector potentials of two electrons. 
In the diagrammatic language of the previous section these latter terms have not been 
considered. The electromagnetic field fluctuations may be decomposed in the longitudinal 
and transverse parts. Here we are interested in the effect of Coulomb interaction and we 
then consider the longitudinal part only. Using the fluctuation dissipation theorem one finds 
the electric field fluctuations for small frequencies as (u < T) 



:2T, 



^ = -q l q^—lmV R (q,uj). 



OJ 



Inserting the retarded interaction from eq. (|T5]) we arrive at 

T 



q 2 e 2 N Dcj 2 ' 



(54) 



(55) 



Note that the fluctuation dissipation theorem has to be used in the classical limit, uo < T. 
As it is clear from the discussion in the beginning of this section it would be inconsistent to 
take the quantum limit 



to 



(A'JI?)^ = coth ) Im^(q, uj) 



(56) 



since this would immediately introduce contributions from uj > T. For example in equation 
(flOl) for the inelastic scattering time this would amount to keep the coth(u;/2T) while ne- 
glecting the tanh(u;/2T). Various claims about the saturation of the dephasing time appear 
to be related to this mistake [j3l,[2~l~|. 

Next we will integrate eq.(B3) by parts by using the relations 



r\r\^q l q 3 exp (iq • (r h -r ta )) = d tl d t2 ex P (iq ■ ( r h - r t2 )) 



(57) 



d tl d t2 fcos(cj^) + cos(c;^ 



cos w 



2 



COS (UJ 



2 

ti-t 2 



(58) 



The integration over the frequency uj may be extended to infinity when considering the 
asymptotic behavior rj ^> 1/T (which involves times t\ and t?. which are much longer than 
1/T) and one finds (^-functions 



1 '" x 
2 



duj h±t 2 , , 
— cos(w — - — ) = 5{ti ± t 2 ) 

Z7T Z 



(59) 



16 



As a result of inserting eqs.([55|-^9|) in eq. ( |53"D one obtains 

T 



Si 



int 



Eq- (pUD is the first central result of this section. It shows that the small q singularity due 
to long range Coulomb interaction is cut off. There is therefore no need to introduce the 
refined form eq. (|4Tf) for the Coulomb interaction. This happens because the oscillating 
term, originating from correlation between different particles, cancels with the unity factor 
stemming from correlations on the same particle. In the diagrammatic language this is the 
cancellation between self-energy and vertex corrections in the small q limit. To proceed 
further, one has to average the oscillating factor over the diffusive paths. However to obtain 
the leading behavior is enough to proceed as follows. The oscillating factor is different from 
zero only when |q • (r(t) — r(— i))| < 1. Then one may substitute the expression in brackets 
with unity with the condition that Dq 2 \t\ > 1 or Dq 2 > r^ 1 . We note that the integral over 
q is divergent in the ultraviolet limit. This is because we took the integral over uj from minus 
infinity to plus infinity. We may correct this error remembering that for typical frequencies 
Dq 2 ~ uj < T and then use this condition as an ultraviolet cut-off in the q integral of eq. (|60|) . 
Then 

T ^ 1 r<t> T 
1 = Y -z I dtl = r d \yiTts (61) 



which yields the AAK result in eq.( 

The dephasing time in the particle-hole channel can be obtained by an analogous calcu- 
lation. The particle-hole propagator is given by 

D? t ,(r, r') = - Dr tl e-( s » +Sl ', (62) 



T 



IV =r' 



with 

5 = //,dt 1 § (63) 
S 1 = -ieftdhh, ■ (A(r tl ,t! + V2) - A(r tl ,*! -tj/2)) (64) 

This differs from the expression for the cooperon by a factor two in So, and what is most 
important, the external field couples differently in Si. 

Averaging over the fluctuating vector potential and going through the same steps as 
before for the particle-particle propagator we find 



^ = \ f dti f dh ? oh~ 2 exp [iq ■ ^ - rtJ] 



x {28{t x - t 2 ) - 5{t x - 1 2 - 77) - 5(h -t 2 + 7])) . (65) 

We set the integration limits as needed in the calculation of Aj EEI , eq. flTTj) , and considered 
the limit r]T ^> 1. Details of the explicit evaluation of the dephasing time in the particle- 
hole channel are reported in Appendix A. The first delta function in (EE) corresponds to 
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correlations of the field for the same particle (or hole), whereas the two other delta functions 
correspond to correlations between a particle and a hole. From the integration limits it 
follows that these do not contribute to Si nt . The low momentum singularity is only cut-off 
if we calculate the field fluctuations more accurately, 

1 K DKq (66) 



Dq 2 qu 2 + {DKq)' 

according to the expression for the screened Coulomb interaction in our calculation of 1/ r inc i 
in eq. (fUD . We find that the dephasing rate in the particle- hole channel equals the 
inelastic rate l/r ine i, and is thus definitely not the same as the dephasing rate in the particle- 
particle channel. 

Our result that the inelastic rate is equal to the dephasing rate in the particle-hole chan- 
nel is non-generic and is due to the dephasing mechanism considered here. The crucial 
approximation involved is that the interaction is local in time. For an interaction with re- 
tardation effects the inelastic scattering rate and the dephasing rate will not be the same. 
With this regard we want to mention that the processes considered here are not the only 
intrinsic processes leading to inelastic scattering and dephasing. Castellani et al |TJj] calcu- 



lated the dephasing rate adding to the diagrams shown in Fig.^| the corresponding Hartree 
diagrams. The latter correspond to exchange of spin-fluctuations. In two dimensions, a con- 
tribution to l/r^ which is linear in T has been found, which corrects the result given here 
by a numerical factor. Beyond electron-electron scattering also electron-phonon scattering 
is a further source of inelastic scattering and dephasing. Often the latter is dominant at 
high temperature with l/r^ oc T 3 . In general any type of low lying excitation with u < T 
that couples to the conduction electrons causes inelastic scattering and dephasing. In the 
zero temperature limit, however, the number of such excitations goes to zero and so does 
the dephasing rate. 



VI. SUMMARY AND DISCUSSION 

We have studied non-linear effects and dephasing of the quantum corrections to the con- 
ductivity. We worked within a time domain representation. In this representation the weak 
localization correction to the current at time t is related to processes where a particle and a 
hole start diffusing at the same time (t — rj) around a closed loop but in opposite directions. 
We succeeded in deriving a new expression for the interaction correction to the current which 
allowed us to calculate the latter in an arbitrary time dependent electromagnetic field. In 
contrast to weak localization, the size of the interaction correction is not related to diffusion 
around a closed loop. In this case the correction to the conductivity can be written in terms 
of the amplitude describing the propagation of a particle and an hole with the particle mov- 
ing from Xj at (t — 2rj) to xj at (t — rj), and the hole propagating along the same path from 
Xj at (t — rj) to xj at t. 

Due to these different processes contributing to WL and EEI we find substantially dif- 
ferent non-linear effects. In a static electric field, no non-linear effects were found for WL 
whereas in case of EEI non-linear effects arise, when the voltage drop over a thermal length 
is of the order of kgT. Note that we did not take into account heating of the sample. Ex- 
perimental evidence for non-linear effects on EEI which are not related to heating has been 
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found several years ago by Bergmann et al |23[ in gold films. Qualitative arguments why 
this effects should exist were given in that work, but a quantitative theory was not worked 
out. 

near 



Recently non-linear effects have been observed in two dimensional systems |24 , 25 



the metal-insulator transition. The physical origin of the metal-insulator transition and the 
metallic phase in these systems is not clear, but one candidate is quantum interference in 
the interaction triplet channel |JT|. In this case one would expect to observe in the metallic 
phase the non-linear effects which we have described in this paper. The sign of the effect 
would depend on the relative importance of the singlet and triplet contributions to the EEI 
corrections P2| . However, from the available experimental data deep in the metallic phase 
we are not able to separate non-linear effects due to heating from those due to interference, 
so at present we cannot provide a quantitative comparison of theory and experiment. 

Besides the effect of a static electric field we examined the quantum corrections in the 
presence of a microwave field. We confirmed earlier calculations on the effect of a microwave 
field on WL and worked out the effects on EEI. A microwave field of arbitrarily low fre- 
quency can destroy EEI efficiently, if the voltage drop on a thermal length is of order of the 
temperature, similar to the non-linear effects in a static field. A high frequency microwave 
field (hu > kgT) is less effective in destroying EEI, i.e. EEI increases when increasing the 
microwave frequency. This behavior has been observed experimentally by Liu and Giordano 
26| in highly disordered gold films. In those experiments the change of the DC conductivity 



with the amplitude of the microwave field was measured, cx(.E'ac) oc AE\ C) and the factor A 
was found to be frequency dependent. In Ref. |26| a frequency dependence of the interaction 
correction to the conductivity was suggested as a possible origin of the experimental results. 
Whereas experimental results are qualitatively in agreement with our theory (see Fig.[|) 
the experimentally observed effect is much larger (by a factor of about 10 3 ) than predicted 
theoretically. 

Vitkalov et al |27|] studied the conductivity of Si-MOS transistors in the presence of a 
microwave field. For electron densities above n s > 10 12 /cm 2 the experimental results agreed 
with the theoretical considerations on weak localization in the presence of a microwave 
field. At lower density, 2 • 10 11 /cm 2 < n s < 10 12 /cm 2 , they found non-heating effects 



which were incompatible with weak localization. From the parameters given in |27] we 
conclude that non-linear effects related to electron-electron interaction become important at 
low density as we demonstrate here below. In a weak microwave field of low frequency u, the 
weak localization correction is quadratic in the electric field and given by [|TJ owiX-Eac) — 
°wl(0) ~ 0.13 • (e 2 / it h)D(eE AC ) 2 uj 2 T^. The interaction correction on the other hand is 
o'eeK-E'ac) — o"eei(0) ~ 0.08 • (e 2 /irh)h 3 D(eE AC ) 2 /(k B T) 3 . Inserting the parameters given 
in Ref. |27j at density n s = 2 • 10 12 /cm 2 , T = 4.2K, = 0.5 • 10~ n s which corresponds to 
H/Ttf, ~ 1.5K, and hu « 0.5K, it is seen that WL and EEI are comparable in size. Note 
the strong dependence of ctwl on the dephasing time, so varying by two varies the cxwl 
by almost two orders of magnitude. EEI dominates if the dephasing time becomes shorter, 
as is the case for lower electron density. In that regime non-linear effects are only weakly 
magnetic field dependent (due to Zeeman splitting in the spin-triplet channel which we 
did not consider in this paper), and proportional to 1/T 3 , so by reducing the temperature 
from 4.2K to 1.7K the non-linearities are expected to increase by more than one order of 
magnitude. Both statements agree with the experimental observations. 
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Finally, we considered dephasing due to intrinsic electric field fluctuations. For WL we 
confirmed the AAK result, leading to l/r^ oc T in two dimensions. In case of EEI we 



found, as it was suggested earlier by Castellani et al |TJ[], that there is dephasing also in 
the particle-hole channel. However in contrast to this earlier suggestion we found different 
dephasing rates in the particle-hole and particle-particle channels, both of which however 
vanish as T — > 0. Mohanty et al |21| observed recently a saturation of the weak localization 
dephasing time in gold wires. It was shown in the experiments that the saturation was not 
due to heating or magnetic impurities. The dephasing rate saturated below T ~ IK to 
values of the order Ujr^ ~ 1-lOmK. We found that the interaction correction is affected by 
dephasing, so when weak localization saturates, also the interaction correction is expected 
to saturate at still lower temperature. In the samples with the highest saturation values of 



I/t^ a saturation of the interaction correction to the conductivity was indeed observed pi 



below T ~ lOOmK. There remains the task to determine the mechanism causing saturation. 
Clearly it cannot be equilibrium electric field fluctuations. Also experimental evidences 
are against the idea of such an intrinsic decoherence mechanism ||33|| . Altshuler et al 



suggested an external (non-equilibrium) microwave field as the origin of the saturation of 
the weak localization dephasing time. The weak localization dephasing time, = tac, and 
the saturation temperature of EEI, T sat , are functions of the amplitude and frequency of 
the microwave field. It was pointed out by Altshuler et al that a microwave field is most 
efficient in destroying weak localization, when the frequency is close its "optimum" value 
i.e., the inverse dephasing time, lo ~ 1/tac, which leads to l/r AC ~ [D(eE AC /K) 2 } l / :i . In 
this case we expect saturation of EEI at kBT sat ~ ^/tao The experimental saturation 
temperature however is larger and roughly found as &BT sat ~ IO/i/Vac- This could be 
explained with a microwave frequency which is not "optimum" . By decreasing the frequency, 
the WL dephasing time becomes larger according to l/r AC fa [D^eExc/'h) 2 '^ 2 ] 1 ^ , whereas 
the saturation temperature of EEI does not depend on frequency. From fcsT sat ~ 10/i/tac 
we estimate the frequency which is consistent with the experimentally observed saturation 
of WL and EEI as hu ~ fc B T sat /10 5 / 2 . 

On the other hand saturation of the resistance could also occur due to heating. In Ref. 
H strong heating is assumed to set in when the voltage drop over the length of the wire is of 
order of the temperature. From our considerations we found a saturation of the resistance 
when the voltage drop over a thermal length is of the order of the temperature. In the 
experiments the thermal length was much shorter than the system size, so the saturation 
of the resistance is most probably due to heating if the conditions of Ref. H are satisfied. 
This however should be checked experimentally. 

The observed saturation of dephasing an resistance could also be caused by some in- 
trinsic processes, other than the equilibrium electric field fluctuations. Recently it has been 
suggested that two-level-systems may lead to a temperature independent weak localization 
dephasing rate |3^J3^] . To assess the agreement with the experiment one has to calculate the 
dephasing rate in the particle-hole channel using for example the formalism we developed 
in this article. This however is beyond the scope of the present paper. 



20 



ACKNOWLEDGMENTS 



This work was supported by MURST under contract no. 9702265437 (R.R.) and the 
TMR program by the European Union (P.S.). C.C. and R.R. acknowledge useful con- 
versation with M. Sarachik and S. Vitkalov. We thank M.Leadbeater for his help in the 
preparation of the final version of the paper. 



APPENDIX A: DEPHASING IN THE PARTICLE-HOLE CHANNEL 

In this appendix we report the detailed calculation of the correction to the particle-hole 
propagator caused by electric field fluctuations. We start from the path integral represen- 
tation of the diffuson, after averaging over the fluctuating field and after the approximation 

D t- V /2,t-3 V /2( r i r ) = ~j £>r t exp(-S - (Sint)d), (Al) 

where So has been defined in section V. The relevant term for dephasing is S int and is given 
by 

H rv rT T D 2 K 2 

Siat -Jo d H J-T 2tt q DN ijj 2 + (Duq) 2 

x exp[iq(r tl - r i2 ) - iuj(fa - fa)] [1 - cos^)] . (A2) 

The cos(u;?y) term is due to the vertex corrections. The expression we gave in eq.(|5"oD is 
recovered when neglecting the frequency dependence of the field fluctuations and considering 
large times r\ only, such that the ^-integration leads to delta functions in time. Here we do 
not rely on these approximations. First we average (|A2j) over diffusive paths, (exp[iq(r tl — 
r t2 )])d k, exp(— Dq 2 \t\ — fa\) and switch to "center of mass" and "relative" variables in the 
time integrations, t + = (ti + fa)/ 2 and t- — fa — fa, leading to 

(S mt ) d -2j o dt + dt. J_ t — ^ — - - {DKq)2 

x exp(Dq 2 \t_ \ - iutJ) [1 - cos(c^)] . (A3) 

Note that the time integration in this expression covers only half of the region of the original 
integral. We cured this by multiplication by two. Then it is straightforward to integrate 
over t_ and t + , 

rri/2 p2t+ 

/ dt+ dt_exp(-/V|i-l 

JO J-2t + 



Re 



V 



iu + Dq 2 \iu + Dq 2 

For large momenta {Dq 2 > 1/rj) the second term is negligible, so (A4|) is proportional to the 



inverse diffusion propagator. The second term becomes important in the low momentum 
region, and cancels the low-g low-u; singularity of the propagator. Explicitely we find 
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f 7]Dq 2 /[uJ 2 + (Dq 2 ) 2 } for Dq 2 > l/r) 

(A4) = < 7][Dq 2 + (1 - cos cor])/ 7]} /uj 2 for TJg 2 < l/r),UJT) > 1 . (A5) 
[ ry 2 /2 for TJg 2 < l/r],ujr) < 1 

Calculating (S^d for large times rj ^> 1/T we find 

(^nt)^^-^^^ 2 ), (A6) 

from the region of the momentum integration with Dq 2 < T and Dhiq > 1/rj. Setting 
(S'mt)d = 1 at r] = 70 reproduces the FA result for the inelastic scattering time, see eq.fflll). In 
the relevant region of the momentum integration (Duq > 1/rj) the electric field fluctuations 
are nearly local in time. As a consequence vertex corrections are not important in this 
region. 

The opposite limit, r] <C 1/T, which is relevant for <C 1/T, may be more interesting 
since only small times 7] contribute to the interaction correction to the current. In this limit 
the vertex corrections are important, and the dephasing time is strongly reduced below the 
inelastic time, since 1 — cos(u;?7) w (uji]) 2 /2, compare (|A2|). However also in this region it 
is essential to take the refined form ([|1]) of the Coulomb interaction in order to determine 
the electric field fluctuations. The vertex corrections do not cancel the low-g singularity of 
ImV^(g). Going explicitly through the algebra again, we find 

T 4 77 4 Dk 2 

{S ^ )d= 2^WN ln T^- (A7) 



The relevant region in momentum space is Dq 2 < 1/rj, Duq > T. From (|A7|) the dephasing 
time in the limit where <C 1/T is determined as 



(A8) 



1 , Dk 2 Y^ / 1 , Dk 2 , a 2nT _ 1/4 
= =T In (24tt 2 N D 1/4 
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FIGURES 
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FIG. 1. Graphical definition of the cooperon (particle-particle propagator) and diffuson (parti- 
cle-hole propag ator). G A andG R are the advanced and retarded Green's functions which in general 
depend on the external electromagnetic field. 
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FIG. 2. Diagrammatic structure of the correction to the Keldysh component of the Green's 
function before (a) and after (b) averaging over impurities. The vertices Ai and A2 are singular 
and depend on the external electromagnetic field. Note the different structure of Ai and A2: 
Whereas the bare interaction vertex in Ai is connected to a retarded and a Keldysh function, the 
bare vertex in A2 is directly connected to a retarded and an advanced function. 
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FIG. 3. Temperature dependence of the conductivity Aceei = ^Jeei/E due to the interplay 
of disorder and interaction in two dimensions. The full line represents the current in presence of a 
static electric field E. The dashed lines represents the conductivity, in presence of a low frequency 
AC field of strength E\q = 10-E which superposes the static field E. The AC field may be parallel 
or perpendicular to the DC electric field. The temperature is in units of T , with D(eE) 2 = T 3 . 




FIG. 4. Temperature dependence of the additional current due to the interplay of disorder and 
interaction (d = 1). The full line represents the current in presence of a static electric field E. The 
dashed line represents the current, in presence of a low frequency AC field of strength Eac = 10.E 
which superposes the DC field E. The AC field is assumed to be parallel to the DC field. We use 
arbitrary units for the current, and the temperature is in units of To, with D(eE) 2 = T$. At high 
temperature, the 1/a/T behavior which is characteristic in one dimension is seen. In this limit, the 
current is linear in the electric field. Below irT ~ To, the current saturates, due to the non-linear 
effects. In presence of the AC field the current saturates at higher temperature. 
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FIG. 5. Conductivity in two dimensions in presence of both a DC field E and a high frequency 
field -Eac of frequency u>. The temperature which is kept fixed is T = 10To, where Tg = D(eE) 2 . 
In the absence of the AC field one is in the linear response regime, according to Fig.^]. The AC 
field supresses the quantum correction to the conductivity, but the effects becomes weaker when 
uj > T. 



(a) 



(b) 



(c) 



(d) 




FIG. 6. Self-energy diagrams for the cooperon and diffuson. The triangles denote impurity 
dressed vertices as in Fig.||. The dashed line with a cross is a single impurity line. Two dashed 
lines denote an impurity ladder, i.e. a cooperon or diffuson. Diagram (e) is not included in the 
calculation of the inelastic scattering, but is included in the path integral approach to the dephasing 
time. 
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